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The possibility of synthesizing a doubly magic superheavy nucleus, ■^^**114i84, is investigated on 
the basis of fluctuation-dissipation dynamics. In order to synthesize this nucleus, we must generate 
more neutron-rich compound nuclei because of the neutron emissions from excited compound nuclei. 
The compound nucleus ■^"^114 has two advantages to achieving a high survival probability. First, 
because of small neutron separation energy and rapid cooling, the shell correction energy recovers 
quickly. Secondly, owing to neutron emissions, the neutron number of the nucleus approaches that 
of the double closed shell and the nucleus obtains a large fission barrier. Because of these two effects, 
the survival probability of '^'^'*114 does not decrease until the excitation energy E* = 50 MeV. These 
properties lead to a rather high evaporation reside cross section. 



PACS numbers: 24.60.Ky, 25.70.Jj, 27.90.-fb 



The search for new elements is a long-standing im- 
portant subject in nuclear physics 0, J2|- According 
to macroscopic- microscopic calculations |3| , there should 
be a magic island of stability surrounding the doubly 
magic superheavy nucleus containing 114 protons and 
184 neutrons. Attempts to synthesize heavy elements 
with atomic numbers beyond Z ~ 100 have been active 
since the 1970s, making use of various developments in 
experimental techniques a i. For superheavy ele- 
ments around Z ^ 114 and N ~ 184, practical combi- 
nations of a target and projectile, such as ^^Ca+^^^'Pu, 
have been used by the FLNR group In this case, the 
neutron number of the compound nucleus is less than 
N = 184. 

Actually, if we plan to synthesize the doubly magic 
superheavy nucleus ^^*114i84, we must fabricate more 
neutron-rich compound nuclei because of the neutron 
emissions from excited compound nuclei. Since combi- 
nations of stable nuclei do not provide such neutron- 
rich nuclei, the reaction mechanism for nuclei with Z = 
114, A'' > 184 has rarely been investigated until now. 
However, because of the characteristic properties of these 
nuclei, we find an unexpected reaction mechanism for en- 
hancing the evaporation residue cross section. We report 
this mechanism in this paper. 

As is well-known, in heavy systems around Z 80, the 
trajectory calculations with friction 0,0 were very use- 
ful for the explanation of the extra- or extra-extra-push 
energy. In superheavy mass region, however, the mean 
trajectory calculations are not suitable, because mean 
trajectories cannot reach the spherical shape region and 
around due to the strong dissipation However, the 
extremely small part of distribution can be found there 
due to fluctuation. Therefore, it is important to take into 
account the fluctuating part from the mean trajectory. 
It becomes necessary to solve a full dissipative dynamics, 
or a fluctuation-dissipation dynamics with the Kramers 
(Fokker- Planck) equation or with the Langevin equation 
llllllll. 

Using the same procedure as described in reference 
|l2j | , we apply the fluctuation-dissipation model and em- 



ploy the Langevin equation for the fusion process. On 
the basis of our previous studies 0, Q , to investigate 
the fission process, we employ the Smoluchowski equation 
which is a strong friction limit of Fokker-Planck equation. 
Here, we take into account the temperature-dependent 
shell correction energy. 

The evaporation residue cross section aER is estimated 

as 



■ ^(2/+l)r,(£;c™, l)PcNiE*,l)WiE*,l), 



1=0 



(1) 

where /io denotes the reduced mass in the entrance chan- 
nel. Ecm and E* denote the incident energy in the centcr- 
of-mass frame and the excitation energy of the compound 
nucleus, respectively. E* is given as E* — Ecm — Q with 
Q denoting the Q— value of the reaction. Ti{Ecm,l) is 
the capture probability of the Ith partial wave, which is 
calculated with the empirical coupled channel model |0| . 
Pcn{E* , I) is the probability of forming a compound nu- 
cleus in competition with quasi-fission events. W{E* ,1) 
denotes the survival probability of compound nuclei dur- 
ing de-exciting process. 

To calculate Pcn, we employ the Langevin equa- 
tion. We adopt the three-dimensional nuclear deforma- 
tion space given by two-center parameterization [l^ Il7| | . 
The three collective parameters involved in the Langevin 
equation are as follows: zp (distance between two poten- 
tial centers), S (deformation of fragments) and a (mass 
asymmetry of the colliding nuclei); a — {Ai — A2)/{Ai + 
A2), where Ai and A2 denote the mass numbers of the 
target and projectile, respectively. The detail is ex- 
plained in reference jlJI- 

After the probabilities reaching the spherical shape and 
around, we must treat extremely small probabilities in 
the decay process of the compound nucleus. Therefore, 
we investigate the evolution of the probability distribu- 
tion P{q, I; t) in the collective coordinate space with the 
Smoluchowski equation We employ the one- 

dimensional Smoluchowski equation in the elongation de- 
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gree of freedom zq, which is expressed as fohows; 



(a) 



dr.,,. 1 <9 (dV(q,i,t) , , T 92 

(2) 

q denotes the coordinate specified by zq. V{q,l]t) is the 
potential energy, and the angular momentum of the sys- 
tem is expressed by Z. fi and /3 are the inertia mass and 
the reduced friction, respectively. For these quantities 
we use the same values as in references [ISIli. T is the 
temperature of the compound nucleus calculated from 
the excitation energy as E* — aT^ with a denoting the 
level density parameter of Toke and Swiatecki . The 
temperature dependent shell correction energy is added 
to the macroscopic potential energy. 



V{q,l;t) = VuM{q) + 



2/(9) 



v,x,Mm. (3) 



where I{q) is the moment of inertia of rigid body at coor- 
dinate q. Vdm is the potential energy of the finite range 
droplet model and ^hcii is the shell correction energy at 

r = o 

The temperature dependence of the shell correction en- 
ergy is extracted from the free energy calculated with sin- 
gle particle energies 14,19]. The temperature-dependent 
factor <i>(t) in Eq. Q is parameterized as; 



^{t) exp - 



(4) 



following the work by Ignatyuk et aZ.|2^. The shell- 
damping energy Ed is chosen as 20 MeV. The cooling 
curve T(t) is calculated by the statistical model code 
SIMDEC 0, ^1 . We assume that the particle emissions 
in the composite system are limited to neutron evapora- 
tion in the neutron-rich heavy nuclei. When the temper- 
ature decreases as a result of neutron evaporation, the 
potential energy V{q^l] t) changes due to the restoration 
of shell correction energy. 

The survival probability W {EQ,l]t) is defined as the 
probability which is left inside the fission barrier in the 
decay process; 



W{E*,l-t) = 



P{q,l;t)dq. 



(5) 



inside saddle 



Here, -Eq is the initial excitation energy of the compound 
nucleus. 

For the purpose of understanding well the characteris- 
tic enhancement in the excitation function, we first dis- 
cuss the evaporation residue probability of one partial 
wave, i.e., of / = 10, which is one of the dominantly con- 
tributing partial waves [l4|. 

In our previous studies [isl IT^ , we showed the isotope 
dependence of the evaporation reside cross section for 
Z = 114. At that time, we investigated the isotope de- 
pendence with N < 184. However, in order to synthesize 



9 

8 

7 

> 

g. 6 

A 

4 

3 





































A^92 


















A=2 


98 


















A= 


304 






















































1 









160 



170 



180 



190 



200 



210 



(b) 



14 £* = 40 MeV 

















































































A=292 
















A=298 
A=304 



















500 1000 

t (10""' sec) 



FIG. 1; (a) Neutron separation energies averaged over four 
successive neutron emissions (Bn) for the isotopes with Z = 
114 [23. (b) Cooling curves of yl = 292,298 and 304 with 
Z = 114 at the initial excitation energy Eq = 40 MeV, that 
are derived by the statistical code SIMDEC [T l ll9 | . 



the doubly magic nucleus ^^^114184 by hot fusion reac- 
tions, we must fabricate a more neutron-rich compound 
nucleus of > 184 because of the neutron emissions from 
the excited compound nucleus. Taking into account the 
neutron emissions, we investigate the possibility of syn- 
thesizing 2^^114i84. 

The neutron separation energy depends on the neu- 
tron number. Figure 1(a) shows the neutron separation 
energies averaged over four successive neutron emissions 
(Bn) for the isotopes with Z — 114. We use the mass ta- 
bic in reference 21]. (B„)s of A = 292, 298 and 304 are 
6.43, 5.25 and 4.06 MeV, respectively. With increasing 
neutron number of the nucleus, the neutron separation 
energy becomes small. Therefore many neutrons evapo- 
rate easily from the neutron-rich compound nuclei. 

Because of rapid neutron emissions, the cooling speed 
of the compound nucleus is very high. Figure 1 (b) shows 
the cooling curves oi A = 292, 298 and 304 at the initial 
excitation energy Eq = 40 MeV, that were derived using 
the statistical code SIMDEC [lj,[ll]. In the case of A = 
304, the excited compound nucleus cools rapidly and the 
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FIG. 2: (a)_Shell correction energies Vshdi of isotopes with 
Z = 114 |2l| . (b) Time evolution of the neutron number 
for the de-exciting nucleus ^°*114i9o for eight different initial 
excitation energies. 



fission barrier recovers at a low excitation energy. 

Moreover, owing to the neutron emissions, the neu- 
tron number of the de-exciting nucleus with A = 304 
approaches that of a nucleus with the double closed shell 
Z = 114, = 184. Figure 2(a) shows the shell correction 
energies Vshei\ of isotopes with Z = 114 T4hen of the 
A = 304 (TV — 190) nucleus is smaller than that of the 
A = 298 {N — 184) nucleus. However, in the de-exciting 
process of the nucleus with A = 304 (A^ = 190), the 
neutron number approaches = 184 because of neutron 
emission. In Fig. 2(b), the time evolution of the neutron 
number for the compound nucleus '^'^*114i9o is shown for 
eight different initial excitation energies, as calculated 
by SIMDEC [ll[l3. At a high initial excitation energy, 
the neutron number of the compound nucleus quickly 
approaches N ~ 184, which is that of a neutron closed 
shell. This means the rapid appearance of a large fission 
barrier. 

The compound nucleus with '^'^'*114 has two advantages 
to obtaining a high survival probability. First, because 
of small neutron separation energy and rapid cooling, 
the shell correction energy recovers quickly. Secondly, 



FIG. 3: Time evolution of the fission barrier height Bf for 
the de-exciting nuclei (a) ^^®114 and (b)^"*'ll4. 



because of neutron emissions, the number of neutrons in 
the nucleus approaches that in the double closed shell, 
and a large shell correction energy is attained. 

Generally, at a high excitation energy, the recovery of 
the shell correction energy is delayed. On the other hand, 
at a low excitation energy, the shell correction energy is 
established. Figure 3(a) shows the time evolution of the 
fission barrier height Bf for ^^^114. We can see that the 
restoration of shell correction energy is increasingly de- 
layed with increasing excitation energy. Using the Smolu- 
chowski equation, we calculate the survival probability, 
which is denoted by the red line in Fig. 4. With increas- 
ing excitation energy, the survival probability decreases 
drastically. 

However, for ^°^114, the situation is opposite. At an 
excitation energy of 50 MeV, the fission barrier recovers 
faster than in the cases with lower excitation energies, 
as shown in Fig. 3(b). The reason is the double effects, 
that is to say, the rapid cooling and rapid approach to 
N 184. The survival probability of ^°''ll4 is denoted 
by the blue line in Fig. 4. It is very interesting that 
the excitation function of the survival probability has a 
flat region around = 20 - 50 MeV. At E* = 50 
MeV, the survival probability of ■^""'114 is three orders 
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FIG. 4: Survival probabilities for ^^s-^^^mo ^-^^ ■^o^lM, 
which are calculated by the one-dimensional Smoluchowski 
equation. 



magnitude larger than that of ^^^114. For reference, the 
survival probability of ■^0''114 is denoted by the green line 
in Fig. 4. 

For the ideal combinations for synthesizing the com- 
pound nucleus ^°'*114, the fusion probabilities for each 
system are shown in Fig. 5, which were calculated 
using the Langevin equation, except for the reaction 
^^^La+^^^La. This symmetric reaction system with ex- 
tremely small fusion probability is applied to the one- 
dimensional Smoluchowski equation. The combinations 
of the projectile and target are indicated in Fig. 5. The 
corresponding Bass potential barriers are denoted by the 
arrows We show the fusion probabilities above the 
barrier, because we use the classical models. To multi- 
ply the fusion probability by the survival probability of 
^°^114 in Fig. 4, we obtain the evaporation residue cross 
section of superheavy elements, as shown in Fig. 6. The 
cross section is rather high. It is expected that neutron- 
rich isotopes are more favorable for the enhancement of 
the evaporation residue cross section. Thus, an investiga- 
tion of the experimental feasibility of obtaining neutron- 
rich superheavy elements is an extremely urgent subject. 

In summary, using the three-dimensional Langevin 
equation for the fusion process and the one-dimensional 
Smoluchowski equation for the survival process on the ba- 
sis of our previous works 0, 0, 0| , we investigated the 
possibility of synthesizing the doubly magic superheavy 
nucleus ^^^114i84. Because of the neutron emissions, we 
must generate more neutron-rich compound nuclei. The 



FIG. 5: Fusion probabilities for each ideal system lead- 
ing the compound nucleus ^'''*114, which are calculated by 
the three-dimensional Langevin equation. For the reaction 
^^^La+^^^La, it is calculated by the one-dimensional Smol- 
couski equation. The arrows denote the corresponding Bass 
potential barriers ,22,|. 
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FIG. 6: Excitation function of the evaporation residue cross 
section for each reaction forming the nuclei with Z — llA. 
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compound nucleus ^"114 has two advantages to achiev- 
ing a high survival probability. First, because of small 
neutron separation energy and rapid cooling, the shell 
correction energy recovers quickly. Secondly, owing to 
neutron emissions, the neutron number of the nucleus ap- 
proaches that of the double closed shell. Because of these 
two effects, the excitation function of the survival prob- 
ability of 114 has a flat region around £;* = 20 ~ 50 
MeV. These properties lead to a rather high evaporation 
reside cross section. As a more realistic model, we plan 
to take into account the emission of the charged particles 
from the compound nucleus. 

Although the combinations of stable nuclei cannot 
yield such neutron-rich nuclei as Z = 114 and TV > 184, 
we hope to make use of secondary beams in the future. 



We believe, the mechanism that we discussed here can 
inspire new experimental studies on the synthesis of su- 
perheavy elements. Also, such a mechanism is very inter- 
esting and can be applied to any system that has the same 
properties, small neutron separation energy and slightly 
larger neutron number than the closed shell. 
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